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ABSTRACT 

We show that magnetized, radiation dominated atmospheres can support steady state patterns of 
density inhomogeneity that enable them to radiate at far above the Eddington limit, without suffering 
mass loss. The inhomogeneities consist of periodic shock fronts bounding narrow, high-density regions, 
interspersed with much broader regions of low density. The flow of radiation avoids the dense regions, 
which are therefore weighed down by gravity, while gas in the low-density regions is slammed upward 
into the shock fronts by radiation force. As the wave pattern moves through the atmosphere, each parcel 
of matter alternately experiences upward and downward forces, which balance on average. We calculate 
the density structure and phase speed of the wave pattern, and relate these to the density contrast and 
the factor by which the net radiation flux exceeds the Eddington limit. The presence of a magnetic held 
is essential for the existence of these flows, since magnetic tension shares the competing forces between 
regions of different densities, preventing the atmosphere from blowing apart. There appears to be a broad 
family of modes propagating in arbitrary directions with respect to the direction of the mean magnetic 
field, and exhibiting a range of density contrasts. While the transition from low to high density occurs 
through a strong shock, the gas must pass through a slow magnetosonic critical point in order to return 
to the low-density state. 

The flux of radiation escaping from the atmosphere exceeds the Eddington limit by a factor of order 
the square-root of the ratio between maximum and minimum density. In principle, this factor can be as 
large as the ratio of magnetic pressure to mean gas pressure. Although the magnetic pressure must be 
large compared to the mean gas pressure in order to support a large density contrast, it need not be large 
compared to the radiation pressure. These highly inhomogeneous flows could represent the nonlinear 
development of the "photon bubble" instability discovered by Gammie. If they occur in nature, these 
structures could have an impact on our understanding of luminous systems such as accreting compact 
objects and very massive stars. 

Subject headings: accretion: accretion disks - black hole physics - hydrodynamics 



1. INTRODUCTION 

The Eddington limit is of fundamental importance to 
the study of luminous systems such as accreting compact 
objects, novae and gamma-ray bursts, and supermassive 
stars. Hydrostatic atmospheres dominated by radiation 
pressure radiate at close to the Eddington limit for the 
central mass, provided that they do not develop density 
inhomogeneities on scales much smaller than the radiation- 
pressure scale height. Accreting black holes either swallow 
all the energy produced in excess of the Eddington limit 
(Begelman 1979), or else the pressure of escaping radia- 
tion constrains the accretion rate to an Eddington-limited 
value (Shakura & Sunyaev 1973). Including the effects 
of rotation leads at best to a modest enhancement in the 
radiation flux (Paczyhski & Wiita 1980). It is possible 
to exceed the Eddington limit by advecting radiation out- 
ward in a very optically thick wind, but this requires a 
kinetic energy flux that greatly exceeds the radiation flux 
(M. Rees 1976, private communication; Meier 1982; Becker 
& Begelman 1986). This is why models for gamma-ray 
bursts demand the re-conversion of kinetic energy back to 
radiative form via external or internal shocks, rather than 
relying on direct emission from the fireball itself (Rees & 
Meszaros 1992, 1994; Narayan, Paczyhski, & Piran 1992; 



Katz 1994; Sari & Piran 1995). 

It has long been known that Eddington's constraint on 
the maximum luminosity of flows and atmospheres can be 
circumvented if the density structure becomes highly in- 
homogeneous on small scales. While diffusive effects tend 
to smooth out small-scale radiation pressure gradients, the 
same cannot be guaranteed of the density if the gas pres- 
sure is comparatively small. The radiation flux, which 
determines the force, is then inversely proportional to the 
local density. Radiation flows readily through tenuous re- 
gions, while avoiding regions of high density. Depending 
on the distribution of high- and low-density regions, the to- 
tal flux through the system can exceed the Eddington limit 
by a large factor (Shaviv 1998, 2000). The matter in the 
low-density regions would be subject to a super-Eddington 
flux, and would be blown upward. But the flux passing 
through the dense regions could be sub-Eddington, result- 
ing in a net downward force due to gravity. At worst, only 
the low-density component of the atmosphere — which 
need not contain most of the mass — would be blown 
away. At best, even the low density component could be 
reined in, either by the inertia of the surrounding dense 
matter with which it interacts, or by magnetic fields cou- 
pling the high- and low-density regions. 

It has remained an open question whether low-density 
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conduits for radiation would form spontaneously in a 
radiation-dominated atmosphere and, if they do, whether 
the coupling between low- and high-density regions would 
be sufficiently strong to allow the atmosphere to main- 
tain overall dynamical equilibrium. Prcndergast & Spiegel 
(1973) proposed an analogy between a radiation-supported 
slab and laboratory "fluidized beds" consisting of com- 
pressed air or water forced through a bed of particles 
like sand. By analogy with the instabilities that occur in 
such flows, they (and Spiegel 1977) proposed that super- 
Eddington fluxes of radiation could escape from radiation- 
dominated slabs via "photon bubbles," i.e., pockets of 
high-entropy (low-density) fluid percolating through the 
higher density background. Recent calculations have pro- 
duced evidence that overstable acoustic modes may indeed 
exist in radiation-dominated atmospheres, although their 
growth rates may be rather small and their robustness un- 
certain (Spiegel & Tao 1999; Shaviv 1999). 

Far more robust unstable modes appear to be present 
when the fluid is permeated by a moderately strong mag- 
netic field. The first photon bubble modes discovered 
(Arons 1992; Hsu, Arons, & Klein 1997) depend on small 
phase shifts due to diffusion effects and are mildly over- 
stable, although they appear to be important in regions 
with extremely high optical depths and strong magnetic 
fields, such as neutron star accretion columns (Klein et 
al. 1996a, b; Jernigan, Klein, & Arons 2000). Gammie 
(1998) later found strongly unstable modes which appear 
to occur under wide-ranging conditions where radiation 
pressure exceeds gas pressure, such as the inner regions 
of accretion disks and the outer envelopes of supermassive 
stars. Both Arons's and Gammie's instabilities require a 
magnetic field strong enough to constrain the fluid motion 
to nearly one dimension. 

In this paper we present a calculation of steady-state, 
inhomogeneous structure in radiation-dominated atmo- 
spheres. Motivated by the linear instability calculations 
of Gammie and Arons, we pay special attention to the 
effects of magnetic fields. Using a local approximation 
confined to scales much smaller than the scale height of 
the atmophere, we derive periodic, plane-wave solutions 
that are fully nonlinear, corresponding to trains of shock 
fronts sliding through the atmosphere. Matter cycles be- 
tween the dense shocked phase, which is pulled downward 
by gravity, and a tenuous phase which is pushed upward 
by radiation force. Magnetic tension coupling the phases 
allows the system to maintain overall dynamical equilib- 
rium. 

The radiation flux passing through the atmosphere is 
inversely proportional to the distance between shocks and 
exceeds the Eddington limit by a factor of order the 
square-root of the ratio between maximum and minimum 
density (~ the ratio between maximum and mean density). 
This factor may be as large as ^> 1, the ratio between 
magnetic and gas pressure. We speculate that these wave 
trains may represent the nonlinear evolution of Gammie's 
(1998) photon bubble instability. 

2. EQUATIONS AND ASSUMPTIONS 

The simplified equations of our model are derived from 
the usual MHD and radiation-hydrodynamic equations, 
the latter assuming the Eddington approximation with 
terms retained only to 0(v/c) in the fluid motions (Mi- 



halas & Mihalas 1984). Instead of modeling the radiation- 
gas thermal coupling (which should have little effect on 
the solutions we seek) in detail, we separate the gas and 
radiation energetics and assume an isothermal equation of 
state for the gas, with sound speed c g = (pg/p) 1 ^ 2 . We 
also assume a uniform gravitational field — gz. The basic 
equations are then 
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where p r is the radiation pressure, F is the radiation flux 
(relative to the fluid frame), and k is the opacity (assumed 
to be constant). By neglecting the entropy of the gas in 
eq. ([}]) we makes a fractional error of order p g /p r , which 
we assume to be negligible to lowest order. We neglect the 
DF/Dt term in eq. (m and simplify eq. (n) by neglecting 
\D\np r / Dt\ compared to \D\n p/Dt\, to obtain 



V • F = -4p r (V • v). 



(7) 



The latter approximation is justified according to eq. (|^) 
if radiation is not strongly "trapped" by the inhomo- 
geneities, i.e., if the condition t> = pn\ < c/t) is satisfied, 
where A is the characteristic separation between density 
enhancements (the wavelength if the flow is periodic) and 
p and v are characteristic densities and velocities, respec- 
tively. In this case the finite rate of radiative diffusion 
produces a force term that behaves like a frictional drag. 

Implicit in the idea of a radiation-supported atmosphere 
is an overall vertical stratification, with a pressure scale 
height |Vlnp| _1 = H ~ p/(pg). To simplify the model 
further, we will seek structures of sufficiently small sep- 
aration that we can neglect this stratification to lowest 
order, i.e., we treat X/H as a first-order quantity. Thus 
we treat p r as a constant to lowest order. The validity 
of all these approximations can be checked a posteriori. 
We then seek "plane wave" solutions in which the vector 
quantities (v, B, F) lie in the x — z plane and all quantities 
depend on position and time through the combination 



s = x + Qz + not. 



(8) 



Surfaces of constant density (i.e., constant s) are tilted 
with respect to the vertical by an angle 6 S = — tan -1 £. 
We also seek solutions that are periodic in s with wave- 
length A. 

Using a prime to denote differentiation with respect to 
s, we write the continuity equation as 



v p' + {pv x y + c(p«*)' = 



(9) 



which integrates to 



vqp + pv x + (pv z = M = const. 



(10) 



By demanding that there be no net mass flux through the 
atmosphere, 



s+A rS-\-X 

pv x ds = / pv z ds — (11) 



and defining the mean density 

rs+X 



pds = p X, 



(12) 



The x— and z— components of the momentum equation 
are, respectively, 

Pov v' x + <? g p' = + ( 2 )A'{b - A) + 



(1 + C 2 )c 



+ (9 )(p-Po) (19) 



p v v' z + Cc 2 gP ' = -^(1 + ( 2 )A'(1 + (A) + 



(1 + C 2 )c 



g)(p-p ). (20) 



we obtain M = povo- The mean density po has special 
significance as the value of the density for which gravity 
and radiation forces exactly balance. It is also the density 
at which the flow must pass through a critical point (asso- 
ciated with the speed of slow magnetosonic waves) while 
moving between regions of high and low density. 

The condition V • B = is automatically satisfied if we 
represent B by 



B X =B(1 + (A); B z = B{b~A) 



(13) 



where B and b are constants and A(s) is any periodic func- 
tion of s with zero mean. Using the continuity equation 
and demanding that A integrate to zero across a wave- 
length, the flux-freezing equation (0) integrates to 



.4 



Povo 



-{bv x 



(14) 



Next we consider the radiation force term in the mo- 
mentum equation. Averaging over a wavelength, the only 
surviving terms in the momentum equation are 



A" 



-pgz+ ^f) ds = -p Q gz+(^F) = 0, (15) 
c / c 



From this point forward, however, it is more convenient 
to work in rotated coordinates oriented parallel and per- 
pendicular to the surfaces of constant s (i.e., surfaces of 
constant density, etc.). We have 



v x + (v z C,v x - v z 

v± = — „.,,„ ; v\t - 



(1 + (2)1/2- 



(1 + C 2 ) 1/2 



(21) 



and 



B(l + bQ . B[C-b+(l + C 2 )A] 
By=- —j, = M | A9N1/9 . (22) 



(i + C 2 ) 1/2 



(i + C 2 )V2 



Note that B± is a constant. We will also find it convenient 
to replace A by the equivalent function 



= c-b+{i + <: 2 )A 

B ± ~ 1 + bC 



(23) 



The rotated components of the momentum equation are 
then 

Povov'u = fi(l + C 2 ) 1/2 0' + (1 + C 2 ) 1/2 .9(P - Po) (24) 



so we write 

F=-^gz + C-%- P o)x + F D . (16) 
up up 

Fjj , which averages to zero and has a divergence that sat- 
isfies eq. (?]), is the portion of the flux associated with 
radiation drag. We obtain the remaining terms in eq. (16) 



by requiring that eq. (15) be satisfied and that the diver- 



gence of these terms vanish. We would like the drag force 
to vanish when v — > and to be linear in v x and v z , thus 
guaranteeing that the average vanishes. If pd{s) is the 
perturbation in the radiation pressure due to drag, then 
eq. (||) gives = (c/ np)p' D (x + (I), while eq. (^) can be 
integrated to give 



4p r 



■(v x + (v z )-(x + (z). (17) 



i + C ; 

Eliminating the velocity dependence via the continuity 
equation, we obtain the following expression for the drag 
force: 

—F D = —— I -{p-p Q )-(x + C,z). 18 
c (1 + C )c 



Po^i + (1 + C 2 ) 1/2 c>' - -f^(l + C 2 ) 1/2 0V + 



ip r KV 



The continuity equation becomes 

(1 + C 2 ) 1/2 pv± = (Po - p)v , 



(26) 



showing that vj_ vanishes where the local density equals 
the mean. 

2.1. Dimensionless Equations 

At this point it is convenient to non-dimensionalize the 
equations. We define: 

^,-^d + C 2 ) 1 / 2 ; »?=A y^ S -§; (27) 
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These scalings are physically motivated. The fiducial ve- 
locity Wq/(1 + C 2 ) 1 ^ 2 ; is t ne pattern speed, and m is the 
associated Mach number for the gas component. The char- 
acteristic length scale is related to the scale height of the 
gas in the gravitational field. The parameter (3 measures 
the ratio of gas pressure to magnetic pressure (smaller than 
the standard "beta-parameter" of plasma physics by a fac- 
tor 2), while the drag parameter £ is the ratio of radiation 
drag force to gravitational force (parallel to the density 
surfaces) for matter moving at the gas sound speed. The 
parallel and perpendicular components of the momentum 
equation become: 



m 2 u'u 



2 1,1 

m + rj 



P 



+ m4(ri- 1), 



(29) 



(30) 



where a prime now denotes differentiation with respect to 
y. We eliminate u± and U|| by using the continuity and 
flux-freezing equations: 



1-r? 
V 

C- 



n 



l + (b 



(31) 



(32) 



Differentiating and substituting into the momentum equa- 
tions, we obtain 



m " 



V 



<jyq' =-()>' 



m 
V 



(33) 



(34) 



Equations (|33|) and (|34|) can be combined to yield an or- 
dinary differential equation for 77' in terms of 77 and <f> plus 
an energy-type integral that gives an algebraic relation- 
ship between 77 and <f>. To obtain the former we eliminate 
(/)', yielding 



V 



1 



/3m/ 



(l-?7) 



1 



(3m z 



m£ 1 



(3m 2 
77 



(35) 



Equation ( p5| ) has the form of a "wind equation" with a 
critical point at 77 = 1 (p = po). In order for a solution to 
pass this point smoothly, the quantity in square brackets 
on the left-hand side must vanish simultaneously. What 
is the physical significance of this critical condition? If we 
reintroduce dimensional variables, the critical condition 
can be written 



= 0. 



1 + C 2 




( 36 ) 

This is just the dispersion relation for magnetosonic waves 
propagating perpendicular to the constant-density sur- 
faces, where vq/(1 + £ 2 ) 1 / 2 is the phase speed. But 



vq/(1 + C 2 ) 1 / 2 is also the speed of matter crossing the crit- 
ical surface in the frame in which the density pattern is 
stationary (since v± = at the critical surface in the "lab" 
frame; see eq. [p6[). Thus the critical point is a true sonic 
point for disturbances mediated by magnetosonic waves. 
Returning to dimensionless variables, the critical condi- 
tion can be written 

2 1 

771 = 

2/3 



1 + 4>l + P ± V(T+ 4>l + P) 2 - 4/3 , (37) 



where <fi CI is the value of <f> at the critical point and the 
± sign corresponds to fast and slow magnetosonic waves, 
respectively. 

To derive the energy integral we eliminate (77 — 1) and 
integrate the resulting expression. The result is 

<f ( Pm 2 \ n ( m 2 \ 

^+m^\l-^—U + P\r ) + — \=£, (38) 

where e is a constant of integration. 

3. SOLUTIONS 

3.1. General Considerations 

We seek flows in which the density varies between high 
values (7/ > 1), where gravity dominates, and low values 
(7/ < 1) where upward radiation force dominates. Since 
we also seek periodic solutions, in the stationary-pattern 
frame the gas must flow from high to low density and back 
again. In order for this to occur smoothly the flow must 
pass through a critical point. It is appropriate to demand 
that the critical condition be satisfied for 77 decreasing from 
high values to low values, corresponding to a smooth tran- 
sition from subsonic to supersonic flow. Since this corre- 
sponds to drj/dy < for vq > 0, we require 



m£ 1 



(3m z 



< 



(39) 



for the case where the magnetic field is "strong" (J3m 2 < 1) 
and the opposite inequality in the limit of a weak magnetic 
field {13m 2 > 1). Since m and £ are positive by assump- 
tion, this imples that (f> must be negative (positive) near 
the critical point, for the strongly (weakly) magnetized 
cases. If the flow is to be periodic in y, there must also be 
supersonic to subsonic transitions. It is easy to show that 
these transitions cannot occur through a smooth crossing 
of the critical point, except possibly for the unphysical case 
m£ = 0, i.e., where there is no radiation drag. According 
to eq. (|37]) , 4> 2 r would have to have the same value at both 
crossings, however, since drj/dy > for the supersonic- to- 
subsonic transition with vq > 0, the sign of the inequality 
( J38| ) would have to be reversed. The only way to accom- 
plish this is to flip the sign of 4> cr (i.e., of By), but this is 
clearlyincompatible with conservation of the energy inte- 
gral (pq) if m£ 7^ 0. Thus, the transition from supersonic 
to subsonic flow generally must occur at a shock disconti- 
nuity. 

The above argument implies that the nonlinear flows 
considered here cannot exist in the absence of a magnetic 
field. In the nonmagnetized limit {f3 — » 00), eq. ( pq ) re- 
duces to 



77 

— (I + 77) = m£, 
77^ 



(40) 
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which does not permit a smooth transition from subsonic 
to supersonic flow (i.e., drj/dy is positive for all y). This 
conclusion is consistent with findings by Gammie (1998) 
and Arons (1992) that magnetic tension is necessary for 
the existence of linear "photon bubble" instability; we 
shall strengthen this analogy later, in § 4.2. Our rejec- 
tion of the nonmagnetized limit also allows us to choose 
a si gn in the expression for the critical Mach number, 
eq. (p7p. Since the regions of high and low density need 
to be coupled dynamically by magnetic tension, we pre- 
sume that the speed in the pattern frame never reaches 
the fast magnetosonic speed and choose the negative root 
for slow magnetosonic waves. Adoption of the minus sign 
in this relation places severe constraints on the value of 
/3m 2 , the parameter that discriminates between "strong" 
and "weak" magnetic fields. First, we note that ftm 2 is a 



2 r at fixed (3. Then, 
min(/3, 1), i.e., the 



monotonically decreasing function of 4 
setting <j) 2 r — 0, we find that (3m 2 < 
field is always "strong" in this sense. 

Now we consider the shock transition. In dimensionless 
variables, the jump conditions across the shock are 



(3 



m 2 A [ - 



(3m A 



Ar] 



A{4> 2 ) 



A0, 







(41) 



(42) 



where A/ = /_)_ — /_ and we denote downstream (up- 
stream) quantities by + (— ) subscripts. We are espe- 
cially interested in the possibility of large density con- 
trasts, rj + 3> 1 3> t]-, in which case the jump conditions 
may be expressed approximately as 



1 - 



ftm z 



V+ 



m 
V- 



1 



1 - 



(3m 2 



(43) 



(44) 



The maximum attainable density contrast 77+ jr\- is lim- 
ited by the strength of the magnetic field, which must be 
able to resist the difference between the upward radiation 
force in the low-density regions and the downward gravi- 
tational force in the regions of high density. Typically the 
field becomes highly distorted when r\- is smaller than 
(3m 2 , or when r/ + exceeds ft^ 1 ■ Thus, we expect that large 
density contrasts will be possible only if ft <C 1, i.e., if 
the magnetic pressure is much larger than the mean gas 
pressure. Note that this does not mean that the magnetic 
pressure has to be of the same order as the radiation pres- 
sure; it may be considerably smaller. In the limit ft <C 1 
the critical condition may be approximated as 



m 2 = 



1 



1 



O(ft). 



(45) 



To construct an explicit solution for the flow for arbi- 
trary ft we solve eq. (35) using eq. ([58]) to determine 4>(r)), 
subject to the critical condition and two additional self- 
consistency conditions. First, our assumption that the 
function A has zero mean over a wavelength A implies 



( — b \ drj 



1 + bCJ rf 



= 0. 



(46) 



Second, the mean density must equal po, implying 



1+ 



(l-n) 



drj 

V 



o. 



(47) 



By satisfying these conditions we determine the value of 
e and guarantee that the shock jump conditions ( f4l|) and 
( |42] ) are satisfied. 

3.2. Strongly Magnetized Case 

3.2.1. Stiff- Wire Approximation 

The method of solution outlined above does not deter- 
mine a unique value for r]^- (or t] + ), but rather allows one 
to calculate r] + in terms of (or vice versa). Thus, for 
specified values of ft, b, and £ it seems possible to find 
a solution for any density contrast up to the maximum 
value beyond which magnetic support fails. In particular, 
for ft <C rj- <C 7T+ <C ft^ 1 , the terms proportional to ft in 
eqs. (|35|) and (Eq) are negligible throughout. The magnetic 
field Imes are nearly straight and <j> « — b)/(l + bQ is 
approximately constant. Because <f> must be negative ac- 
cording to eq. (|3^), b — £ and 1 + 6£ must have the same 
sign. We refer to this limit as the "stiff-wire approxima- 
tion" because the tension in the magnetic field essentially 
reduces the flow to one-dimensional motion. This approx- 
imation is also relevant to the photon bubble instabilities 
described by Arons (1992) and Gammie (1998). 

In the stiff-wire approximation, the critical gas Mach 
number m reduces to 

|! + & CI 



(l + 6 2 ) 1 /2(l + £2)1/2 

and the wind equation assumes the simple form 



if 

-r(l + ?7) =m£ 

7]* 



fr-C 
1 + 6C" 



(48) 



(49) 



The condition ( |39| ) places an upper limit on the drag pa- 
rameter £ consistent with a solution: m£ < (b—()/(l + b(), 
corresponding to 



4p rKCg ^ (i + c 2 )(i + fo 2 ) 1/2 (6-C) 



< 



(50) 



g c (1 + &C)|1 + 1<| 

If the drag is very strong, i.e., if £ 3> 1, solutions are pos- 
sible only for b w — 1/£, i.e., for \Bn\ ^> \B±\. The integral 
condition 



Jjj is readily evaluated, yielding 

V-V+ = 1) (51) 
and the wind equation is easily integrated to yield 

b-C A 11, (r,+ 
— — ~mi\y= h In — 

i + K J v v+ V v 

where y is the scaled distance downstream of the shock. 
The wavelength A (i.e., distance between shocks) is given 
by setting 77 = r;_ = in cq. (j52|) : 



(52) 



A = 



6-C 



1 

V+ 



2 In rj+ 



(53) 



.1 + &C 

Note that, in the submagnetosonic region between the 
shock and the critical surface, the density drops exponen- 
tially with y, as expected for an isothermal, gas pressure- 
dominated atmosphere in a uniform gravitational field. 
This is an important clue to the physical mechanism driv- 
ing the modes. 
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3.2.2. Weak Field Distortion 

At the next level of approximation we can calculate 
small distortions in the magnetic field caused by gravi- 
tational and radiation forces acting on the gas. Adopting 
the ordering 

(3 < ^! < 1, (54) 

n- 

we write — 0o + (50 where 0o — (( — b)/(l + bQ; and 
m = mo + <5m, where m\ = (1 + 0q) _1 - We linearize the 
energy integral (|3^) and the critical condition ( |37|) with 
e = Eq + Se, keeping terms to first order in (3. We then 
apply the integral condition ([Sf ) , using (^) for 77' and the 
unperturbed value of 77+ = I/77-, since the integrand is 
already of first order. In the integral, we keep only the 
leading term in powers of l/?7- ^> 1. By evaluating the 
perturbed critical condition (at 77 = 1, of course) and the 
energy integral at both 77 = 1 and 77 = r/_ , we finally obtain 



5m = —(3 



1 - o m o A O ml 
0o + ™o£ J 2?7_ 



(55) 



m £ 



mg(l 



0oWo£) 



V 



(56) 



1 I 

It is straightforward to determine Se as well. 

Using eq. (^6|) one can verify that the jump condition 
( |43| ) is satisfied to first order. One can use the second 
jump condition, eq. Q), to determine the effect of small 
field distortions on the density contrast. For given ?7_, 
77 + = l/v7_ + <5r7_|_ , where 



11+ 



j / 4 + 0g -30 o m o A 0o m^ 



00 + «i £ 



2m 



3.2.3. Flows with Strong Field Distortion: Limits on 
Density Contrast 

The weak distortion approximation breaks down when 
<50 ~ 0o- For arbitrary values of £, b <~ 0(1) this will oc- 
cur when 77 + /?7_ ~ 0(/3~ 2 ). Since (3m 2 jr\- ~ O(l) in this 
limit, one must formally retain all terms in equations (35) 
and ([38|) and the problem becomes analytically intractable. 
Note, however, that variations of the field with position are 
confined to the region close to the density extremes (i.e., 
close to the shock). In the region (3 <C 77 <C (3~ x is ap- 
proximately constant. To see this we write the solution of 
the quadratic (E|) as 



-m£ 1 



(3m 2 



(58) 



lm 2 e 1 



(3m 2 



2s - 2/3 77 



n 



where the minus sign is dictated by condition (|39|). The 
argument of the square-root must be positive everywhere, 
therefore it is insensitive to 77 except (possibly) near the 
boundaries. 

Although the magnetic field lines are straight (i.e., w 
const.) over much of the flow, we may not use the ap- 
proximation w (£ — b)/(l + bQ when (3m 2 /r]- ~ O(l). 
Instead the integral condition ( plq ) is satisfied by large vari- 
ations in 0(77) close to the limits of integration. To calcu- 
late these variations requires a full numerical treatment, 



which we reserve for a later paper. However, one can place 
upper bounds on the value of /3m 2 /?7_, and thus on the 
maximum density contrast, without solving for the flow 
explicitly. There are at least three separate constraints. 

First, must be real throughout the flow. In order for 
to be real at t?_ <C 77 -C 77+ we must have 2e + m 2 £ 2 > 0. 
The condition for 0_ to be real is then 



(3m 2 



< 



I 



rj- m 2 t; 2 
which imposes a constraint for 

e < 1 



1 + m 2 £ 2 - v / l + 2m 2 £ 2 (l-e)l , (59) 



2m 2 £ 2 



(60) 



Second, the quantity in square brackets on the left-hand 
side of the wind equation ( |35|) must not change sign any- 
where except at the critical point 77 = 1. For 77 <C 1 this 
corresponds to 

^ ' < ■ ^. (61) 



< 1 



i] 



Physically, this amounts to saying that the kinetic energy 
in the flow (which dominates for 77 <C 1) cannot exceed 
the total magnetic energy. Third, a necessary condition 
for the jump condition (J44|) to be satisfied is 



f 



1 



(3m 2 
2^1 



> 0. 



(62) 



This condition essentially guarantees that the component 
of the magnetic field parallel to the shock front is not so 
strong that it quenches the shock. Conditions ( |6l| ) and 
(p2) can be combined into a single condition: 



(3m 2 



77-(l + 0i) 



< mm 



[1,20: 



3.3. Super-Eddington Luminosities 



(63) 



Equation (|52j) indicates that y is inversely proportional 
to 77 for 77 <C 1. Since the vertical radiation pressure 
gradient is uniform to lowest order, the radiation flux 
cx — Vp r / p passing through the gas is also inversely pro- 
portional to the density. Therefore, regions of different 
densities contribute to the total luminosity in proportion 
to y/rj cx ?7~ 2 . The total flux is strongly dominated by the 
regions of lowest density, as pointed out by Shaviv (1998). 

To state this quantitatively, the "mean Eddington fac- 
tor" (the ratio of the spatially averaged flux to the Ed- 
dington flux eg I k) is given by 



dr) 



where the wavelength (in units of y) is given by 

r 1 "- d V 



A 



77' 



(64) 



(65) 



To leading order in the stiff-wire approximation, the ex- 
pression for the Eddington factor is remarkably simple: 



I 



2?7_' 



(66) 
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independent of the orientation of the constant-density sur- 
faces, magnetic field direction, or strength of radiation 
drag. Flows supporting large density contrasts can ex- 
ceed the Eddington limit by a factor of order 1 jr\- 1 
without blowing away. In the strongly magnetized case 
this factor may be as large as ~ the ratio between 

the magnetic pressure and the mean gas pressure. 

4. DISCUSSION AND CONCLUSIONS 

4.1. Physical Nature of the Flows 

To understand the physical effects driving the flow, it 
will suffice to consider the "stiff- wire" regime, i.e., far 
from maximum density contrast. To lowest order, the 
gas slides along straight magnetic field lines in response 
to the force components parallel to the field. In the sub- 
sonic (dense) region between the shock and the critical 
surface those forces are dominated by the actual grav- 
ity, g = —gz, plus a constant "effective gravity" due to 
the density-independent portions of the radiation flux, as 
given by equations (Qq) and (jig): 



goff = 9 [(C + m 0* + Cm^z] 



(67) 



Projecting g+g c ff onto the magnetic field lines and solving 
the equation of hydrostatic equilibrium under the condi- 
tion that p depends on s = x + (z, we obtain 



p oc exp 



(68) 



which agrees exactly with eq. (52) in the limit 1 <C rj < 
Thus the subsonic portion of the flow resembles an isother- 
mal atmosphere supported partially by the magnetic field. 
The density declines exponentially downstream of the 
shock, with an e— folding distance of order the "gravita- 
tional" scale height of the gas. 

The transition to supersonic flow occurs when the ex- 
ponential decline has reduced the density to ~ po, allow- 
ing "radiation force" (those parts of the flux that are in- 
versely proportional to density) to take over from "grav- 
ity". In the supersonic regime, p <C pa, and the den- 
sity decreases more gradually, p oc 1/s, reaching p_ and 
the next shock after traversing approximately 1/??- grav- 
itational scale heights. To conserve mass flux, the speed 
in the supersonic region increases approximately linearly 
with distance. 

To further clarify the nature of the flow pattern, consider 
the simple case in which the density surfaces are vertical 
(C = 0) and radiation drag is negligible (£ = 0). Then the 
subsonic flow is dominated by the real gravitational field 
and condition ( |39| ) implies that b = B z /B x > 0. Since 
we have assumed vq > the pattern of shocks and rar- 
efactions moves toward the left (— x direction), while the 
magnetic field lines slope upward toward the right. Thus, 
the pattern of inhomogeneities runs "downhill" along the 
field lines, rather than uphill, as a result of the combined 
action of radiation pressure and gravity. This makes physi- 
cal sense: the dense regions just downstream of the shocks, 
which are affected more by gravity than by radiation pres- 
sure, are pulled downhill at nearly the pattern speed. Con- 
versely, gas in the low-density upstream regions, affected 
more by radiation force than gravity, are pushed uphill 



into the shocks by the super-Eddington flux. These fea- 
tures of the flow are illustrated schematically in Figure 1. 
The only difference between this simple case and the more 
general situation with nonzero £ and £ is that the "grav- 
ity" force in the latter case includes a contribution from 
the radiation field. When this is taken into account, the 
material behind the shock still runs "downhill" and inter- 
cepts material being blown upward by the radiation force, 
in such a way that upward and downward forces acting on 
any parcel of matter average to zero over time. 

Now let us consider weak distortions of the magnetic 
field. For the rest of this section we consider only the 
strongly magnetized limit. Equation (|5^) combined with 
eq. ( |39| ) shows that 5<j) is positive close to the shock (i.e., 
where 2rj+2m 2 /r] > m 2 jr\-) and negative elsewhere. Since 
<p = B\\/ B±_ and since (f>o is negative, this implies that the 
field orientation becomes more perpendicular to the con- 
stant density surfaces in the vicinity of the shock, and more 
parallel to the density surfaces elsewhere. The physical 
significance of this result becomes clear if we revert once 
again to the simple case £ = £ = 0. In this case 4> = — b 
so a negative value of 5<j) increases the rightward slope of 
the field, while positive 5<p diminishes it. Thus we see that 
the slope decreases on both sides of the shock. This is due 
to the magnetic field "sagging" under the weight of the 
dense, shocked gas. Elsewhere, the field actually steepens 
as the radiation force pushes the gas upwards. The critical 
point lies in the region of steeper field, so it is no surprise 
that the critical Mach number [ss (1 + ^ r ) -1 / 2 ] is smaller 
than it would be in the absence of field-line distortion, i.e., 
Sm < according to eq. ( |55| ) . 

The limit of negligible radiation drag (£ = 0) allows us 
to derive some interesting approximations that are valid 
for strong distortions of the magnetic field. First, we can 
apply the conditions ( |59| ) , (^0|) , and ( |63| ) to obtain bounds 
on the minimum possible density, in the form: 



(3m 2 
V- 



< 



e 

l+2e 



1 + f - iv/e(e-4) 



< e < 1 

1 < £ < 4 
e > 4. 



(69) 



Note that e must be positive in this limit and that a_ 
can never exceed 3, with the maximum possible value oc- 
curring for e = 4. Since m 2 m (1 + 2s)" 1 , the minimum 
possible value of rj- continues to decrease for increasing 
e, but as we show below (eq. fl73]| ), in this limit the value 
of e is determined by the orientation of the magnetic field 
lines and the wavefronts. 

By expressing the (exact) critical condition in the form 



m 2 (l + 2s) = 1 + (3m 2 + 3/3m 4 



(70) 



we can extract a factor (1 — rj) from the left-hand side of 
the wind equation, which can be canceled with the same 
factor on the right-hand side to yield 



(3m 2 



3(3m 4 
V 



2e-2(3[ V + — 



1/2 



(71) 

For large density contrasts the integrals for A, Edding- 
ton factor £, and the normalization condition (|4(]) relating 
e to b are all strongly dominated by the low-density re- 
gions. This is no surprise considering that 1) the dense 
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Fig. 1. — Schematic illustration of the flows discussed in this paper. The shock front makes an angle tan -1 £ with the vertical and slides 
toward the left, intercepting matter blown upward and to the right by radiation force. Between the shock and critical surface matter flows 
downward and toward the left. Depiction of a magnetic field line illustrates the kinds of distortions caused by forces acting on the gas. 



slabs must be much narrower than the interslab regions, 
to conserve mass (since the mean density must equal po); 
and 2) most of the radiation flows through the low-density 
regions. Keeping only leading terms in this regime and 
integrating from r;_ to infinity rather than 77+, we obtain 



A 



1 r 



2e+ (a- - 1)(2e- 2a_) 



1/2 



(72) 



where we have substituted ra 2 m (1 + 2e) 1 where appro- 
priate. Equation (Eq) gives 



a_(l + 2e)- 



U 2 _ 



1 + b( ~ V27 + (a_ - l)(2e - 2a-) 1 / 2 ' 



(73) 



enabling one to relate the relative orientation of the mag- 
netic field and the constant density surfaces to the "energy 
parameter" e. A similar expression can be derived for the 
Eddington factor I. All of these expressions reduce to the 
correct limits in the stiff- wire approximation, a_ « 1. 

4.2. Relationship to Gammie's Photon Bubbles 

There are tantalizing clues that the steady-state wave 
trains we have described may result from the nonlinear 
development and saturation of Gammie's (1998) photon 
bubble instability. The physical mechanisms that sustain 
these nonlinear inhomogeneities, described in § 4.1 above, 
appear to be similar to the physical effects responsible for 
driving Gammie's instability (see especially § 6.3 of Gam- 
mie 1998). In particular, Gammie's instability depends on 
1) approximately one-dimensional motion enforced by the 
tension in the magnetic field; 2) the inverse relationship 
between local radiation flux and density; and 3) propa- 
gation of the pattern of density inhomogeneities downhill 
along field lines. Like Gammie, we have found solutions for 



a wide range of wavelengths, propagation directions, and 
mean magnetic field orientations. Moreover, like Gammie 
we can show that the inhomogeneities go away in the ab- 
sence of a magnetic field. 

Our model does not incorporate the physics necessary 
to describe nonlinear developments of Arons's (1992) in- 
stability. Arons neglects the dynamical effects of gas pres- 
sure, an approximation that amounts to taking the limit 
m 2 — > 00 in eq. (|35|), with j3m 2 and m£ finite. Equations 
( |35| ) and ( pcf ) clearly do not permit solutions that cross the 
critical point in this limit. The development of Arons's in- 
stability is, in fact, driven by the small but finite value of 
Dp r /Dt, which we have neglected. 

4.3. Validity of Approximations 

We made a number of approximations in deriving the 
wave structure. First, we did not consider the thermal 
coupling between the gas and radiation, assuming instead 
an isothermal equation of state for the gas. This approx- 
imation should be rather accurate, given that fractional 
perturbations in the radiation pressure (and hence in the 
LTE radiation temperature) are much smaller than per- 
turbations in the gas density. The waves are sufficiently 
slow that the gas should have time to come into thermal 
equilibrium with the radiation locally. Such thermal ef- 
fects were also found by Gammie (1998) to have negligible 
effect on the linear photon bubble instability. 

Second, we neglected DF/Dt and Dp r /Dt in the radia- 
tion hydrodynamic equations. These assumptions bracket 
the optical depth associated with the wave motion: 



v . c 

- < np\ < -, 

C V 



(74) 



where A represents the scale length of a region with density 
p and v is the associated velocity relative to the pattern 
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frame. The left-hand inequality corresponds to the neglect 
of DF/Dt and is most stringent in the supersonic region 
of the flow, where np\ ~ npoc 2 /g is approximately con- 
stant. Ignoring geometric and numerical factors of order 
unity (and assuming £ ~ O(l)), this condition becomes 
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npoCgC 



(75) 



When condition ( |7q ) is violated the diffusion of radiation 
may be "flux- limited" (Mihalas & Mihalas 1984), although 
the Eddington approximation (i.e., approximate isotropy 
of the radiation field) is still likely to be accurate. Im- 
posing condition ([TS]) on our flows would place a severe 
constraint on the maximum possible value of the Edding- 
ton factor £ (through eq. [p6[) and a strong lower limit on 
the permissible optical depth of the atmosphere (since we 
seek solutions with f)_ < 1). In fact, both of the specific 
cases we consider in § 4.4 may violate this condition. How- 
ever, in practice this limit is irrelevant for the solutions we 
have derived. This is because Dp r / Dt (i.e., radiation trap- 
ping) is negligible when condition ( ]7q ) is violated, hence 
the radiation force term in the momentum equation ^) 
can be calculated without knowledge of the relationship 
between F and Vp r . Our flows are therefore unaffected by 
flux-limited diffusion. 

The right-hand inequality of (|74|), which corresponds to 
neglecting Dp r /Dt and thus radiation trapping, is most 
stringent in the region of exponential density gradient just 
behind the shock front. This condition can be written as 



< 1. 



(76) 



Additionally, our assumption that both p r and Vp r are 
constant to lowest order amounts to sampling a region 
of the atmosphere much smaller than the radiation pres- 
sure scale height, H ~ p r /(pog). This means we can only 
treat wavelengths A -C H. For waves with the maximum 
density contrast, the wavelength is of order (3~ 1 H g , where 
H g = c 2 /g is the scale height associated with the gas. Thus 
our models should be valid for the entire allowed range of 
wavelengths provided that the magnetic pressure is much 
larger than the mean gas pressure but smaller than the 
radiation pressure. For magnetic pressure comparable to 
or greater than the radiation pressure, the global structure 
of the atmosphere would impose a serious bound on the 
maximum wavelengths and density contrasts. 

Finally, Gammie argues that his WKB analysis is 
strictly valid only for wavelengths smaller than HM{j, 
where Mq ~ c/ (upog^^H 3 / 2 ) is a radiative diffusivity pa- 
rameter. In our notation this constraint corresponds to 
A < H g /t; 2 , if we associate the wavelength of a nonlinear 
wavetrain with that of a linear photon bubble mode. We 
have not found any special significance to this scale in the 
nonlinear steady-state analysis. 

4.4. Applicability to Accretion Flows 

In this section we briefly assess the applicability of our 
model to two important accretion problems, taking into 
account the effects of inhomogeneities self-consistently. 

4.4.1. Thin Accretion Disks 



Consider first a thin, Keplerian, radiation-dominated 
accretion disk. For an a-model viscosity, the dissipa- 
tion rate yields a radiation flux F ~ ^ap r Qh, where 

£1 = (GAf/r 3 ) 1 / 2 is the angular velocity and h is the disk 
semi-thickness. Assuming a vertical gravity g ~ Q 2 h and 
a disk surface density E ~ p r /g we obtain the Eddington 
factor 

t~l—JXlh. (77) 
2 c 

Recall that £ is the ratio of the flux to the Eddington 
flux associated with the vertical component of gravity, 
Fe — cg/n. Thus, a value of £ 1 does not necessar- 
ily imply that the disk is radiating a super-Eddington lu- 
minosity for the central mass, but rather that the disk 
is much thinner (by a factor ~ £~ l ) than it would be if it 
were radiating the same flux but without inhomogeneities. 
Assuming LTE, p r = aT 4 /3, we find that the ratio of gas 
pressure (= pkT j \i) to radiation pressure is given by 



Pr 



ON V 4 Z, 



(78) 



If the viscosity is magnetic in origin with magnetic pres- 
sure pm ~ c<Pr, then the magnetic parameter is given by 
j3 ~ Pg/pM ~ ce~ 1 p g /p r . This leads to an important upper 
limit on £, 



£<(3~ 



Pr 
I . 

Pg 



(79) 



The maximum possible wavelength ~ ah is comfortably 
shorter than the scale height. 

Let us now evaluate the various constraints numeri- 
cally for a central mass of toMq, in terms of the radius 
r = xGM/c 2 — 1.5 x 10 5 a;77i. We express E in units of g 
cm j which corresponds roughly to the electron scatter- 
ing optical depth. The upper limit on £ from ( |7S| ) becomes 



< 1.4 x 10 5 aE~ 1/4 



7/4 



-V2 m l/4 



while the expression for energy generation becomes 
£~ 0.6oE (-) 2T 1/2 . 



Eliminating E between (|30|) and (|8l| ) we obtain 
£< 1.2 x 10 4 a( x- 1 / 2 ™ 1 / 5 . 



(80) 



(81) 



(82) 



Note that hjr < 1, implying an upper limit to £. A similar 
argument gives an upper limit to E: 



fh \ 3/5 

E<1.9xl0 4 f-J m 1/5 gcnr 2 . 



(83) 



From the self-consistency condition (|76|) it is straight- 
forward to show that radiation trapping is not important 
for 



E < 8 x 10 



13/11 



(84) 
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By combining conditions ( |83| ) and (|8J) we can show that 
radiation trapping in the inhomogcncities is unimportant 
for 

- > 1.6 x 10-V/ 16 ™- 1 / 8 . (85) 



If ( pa ) is not satisfied, then radiation trapping could be 
important in the densest parts of the flow, necessitating a 
modified treatment. 

Condition (f75|) becomes 



-15/8 



£<5.3x 10- 4 £ 9 / 8 [ ■- ) x^m- 1 '*, (86) 



implying that flux-limited diffusion may occur in the su- 
personic portions of our flows. However, as noted earlier, 
this will have no effect on the dynamics of the solutions or 
the escaping flux. Finally, we note that the radiation drag 
parameter may be expressed as 



£ ~ 4kS ~ 4.3 x 10- 3 £ 9 / 8 (1 



1/8 



x -^m- 1 /* 



4/5 



x-^rn}/^) 



< 280 I - 

r 



where the last inequality comes from condition (|83j) and 
we have neglected the dependence on £. 

4.4.2. Accretion onto Neutron Star Polar Cap 

We next consider supercritical accretion onto the mag- 
netic pole of a 1 Mq neutron star with radius i?* = 10 
km. The surface (Newtonian) gravity is g — 1.3 x 10 14 
cm s -2 and we assume Thomson opacity k = 0.4 cm 2 
g , neglecting possible modifications to the opacity due 
to a strong magnetic field. Energy is deposited at a rate 
£Fe = leg jn = 1.0 x 10 25 1 erg cm -2 s _1 , where Fe is the 
Eddington flux. Note that the total radiatied luminosity 
is then £Le(Q/4:"x), where f2 here is the solid angle sub- 
tended by the accretion mound. To avoid complications 
(such as atmospheric pressure gradients) that are outside 
the scope of our simple model, we consider only the re- 
gion just below the accretion shock, which we assume lies 
<; i?* above the star's surface (Arons 1992). The pressure 

and density are then p r ~ 2£F E /vs{R*) — 1.2 x 10 15 ^ erg 
cm -3 and p ~ 7p r /vg(R*) = 3 x 10~ 5 ^ g cm~ 3 , where 
we have assumed a radiation-dominated adiabatic shock. 
(Note that the pressure and density may be considerably 
higher closer to the star.) The associated temperature 
and gas sound speed in LTE are T = 2.6 x 10 7 i 1 ^ K, 
c g = 6.6 x lO 7 ^ 1 / 8 cm s~ 4 . The ratio of radiation pres- 
sure to gas pressure is lO 4 ^ 1 / 4 , and the drag parameter 
is £ = 0.03^ 9 / 8 (l + C 2 ) -1 / 2 . Radiation drag becomes im- 
portant only if t J> 30, but even then it will not inhibit 
the development of inhomogcncities nor strongly affect 
the maximum luminosity. Although our drag parameter 
£ may be ^> 1, the diffusion parameter of Arons 's anal- 
ysis, Mo ~ (pr/pg)^ 1 / 2 , is <C 1, consistent with his 
approximation scheme. 

For accretion onto a neutron star the magnetic field is 
anchored in the neutron star's crust. We expect the mag- 
netic pressure to exceed the radiation pressure, implying 



that the stiff-wire approximation applies and the maxi- 
mum wavelength of the inhomogeneitics is limited by the 
scale height of the atmosphere, which we take to be h <; 

ii*. This implies (very roughly) I < 2 x 10 3 (/i/10 6 cm) 4 / 5 . 
Trapping of radiation is not an issue here since condition 
@ is satisfied for £ < 2.5 x 10 4 . 

Using eq. (66) to relate ry_ to I, we find that condition 
(|7q ) becomes f 1 / 8 3> 11, which is unlikely to be satisfied. 
However, once again we note that flux-limited diffusion has 
no effect on our flow solutions when radiation trapping is 
unimportant. 

4.5. Summary and Prospects for Future Work 

We have demonstrated that magnetized, radiation- 
dominated atmospheres can support periodic, steady-state 
trains of shocks with wavelengths much shorter than the 
radiation pressure scale height. The most important re- 
sult of our analysis is that the density inhomogeneities 
associated with these shock patterns can allow radiation 
to leak out of the atmosphere at a rate that far exceeds the 
Eddington limit, without blowing the atmosphere apart. 
There are reasons to suspect that these solutions represent 
the nonlinear outcome of Gammie's (1998) photon bubble 
instability. 

The amount by which the radiation flux can exceed the 
Eddington limit is limited by the ratio of the magnetic 
pressure to the gas pressure. It is still an open ques- 
tion whether dynamo action in a radiation-dominated at- 
mosphere can amplify the magnetic pressure to a value 
that greatly exceeds the gas pressure (Sakimoto & Coro- 
niti 1989; Miller & Stone 2000). In neutron star accretion 
columns the magnetic pressure would be fixed by internal 
currents in the star and could greatly exceed the radia- 
tion pressure, let alone the gas pressure. The maximum 
possible radiation flux might then be limited by the scale 
height of the atmosphere. 

We have not yet shown that well-defined, persistent 
wavetrains can exist in real atmospheres. Although our 
solutions are fully nonlinear, they are not global. The so- 
lutions formally admit arbitrary density contrasts (with 
corresponding wavelengths) up to a maximum imposed by 
the magnetic field strength, as well as nearly arbitrary di- 
rections of propagation with respect to the mean magnetic 
field. Our analysis affords no indication of whether there 
are preferred wavevectors or wavelengths, which would be 
essential for establishing persistent structure. If preferred 
scales and directions exist, are they set by the stability 
of the nonlinear solutions (about which we know nothing 
yet), by spatial boundary conditions associated with the fi- 
nite scale height of the atmosphere, or by initial conditions 
associated with the development of these inhomogeneities 
(e.g., the fastest growing modes of Gammie's instability)? 
The entire class of flows may prove to be unstable, as 
one might worry given that they involve regions of low 
density accelerating into high-density slabs. However, the 
fact that the interface is a shock makes it less likely that 
Ray leigh- Taylor instabilities would disrupt these flows. 

We hope to address these issues soon. Ideally, radiation- 
hydrodynamic computer simulations will be able to handle 
these situations within the next year or two (J. Stone 2000, 
private communication); Hsu et al. (1997) have already 
successfully modeled the development of Arons 's photon 
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bubble instability using simplified dynamics. There is also 
further scope for analytic techniques to address stability 
by considering weak time dependence, and the effects of 
boundary conditions by modeling the role of the atmo- 
spheric pressure gradient more accurately. 

These questions seem well worth pursuing. Shock trains 
could be important in systems as diverse as accretion disks 
around black holes and neutron stars, the envelopes of very 
massive stars, and debris from stars that have been tidally 
disrupted by massive black holes in galactic nuclei (Rees 
1988). The generalization of these solutions to radiation- 
dominated outflows could lead to applications in novae, su- 



pcrnovae, and gamma-ray bursts. This phenomenon might 
also permit the rapid growth of supermassive black holes 
by accretion at high redshifts. If radiation-dominated 
atmospheres spontaneously develop such strong inhomo- 
geneities that they can transmit radiation fluxes far in ex- 
cess of the Eddington limit, it would have a qualitative 
impact on our understanding of all such systems. 
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